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ABSTRACT 
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The guiding center motion of particles in a 'nearly' drift free magnetic 
field is analyzed ir. order to investigate the dependence of mean drift velocity 
on equatorial pitch angle, the variation of local drift velocity along the 
trajectory and other properties. Consider a 2-dimensional magnetic field li , 
somewhat resembling the field of the geomagnetic plasma sheet, which has the 
following properties: (l) No y component and no y dependence (2) B x = 0 in the 
equatorial plane z = 0 and (j) B z = B Q = constant. In such a field, , it 
has been shown, particles may be trapped, but they exhibit no net drift in the 
y direction: their instantaneous drifts at any time may be large, but when they 
are averaged over a sufficiently long tine (or, for adiabatic particles, over 
a bounce period) the result is zero. Here a slight modification of this mode is 
explored -- motion in a field in which B slowly varies in the x direction, 
as is indeed observed in the geomagnetic tail. The mean drift <v S for 
adiabatic particles can now be expressed by means of elliptic integrals: it 
no longer vanishes but is merely small -- if v Q is the drift velocity for 
equatorial particles (which Is easily derived), <v S is typically between 
v Q and v q /3 . By contrast, instantaneous guiding center drift velocities at 
z ~ 0 may be 50 times larger. Explicit approximations to the twice -averaged 
Hamiltonian W(o(, p, J) near z = 0 are also derived, permitting simple 
representation of drift paths if an electric potential ) also exists. 

In addition, the use of W or of expressions for the longitudinal invariant 
allows the derivation of the twice averaged Liouville equation and of the 
corresponding Vlasov equation. Bounce times are calculated (using the drift-free 
approximation), as are instantaneous guiding center drift velocities, which are 
then used to provide a numerical check on the formulas for 
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It has been shown £ Stern and Palmadesso , 1975J that In two-dimensional 
magnetic fields of the form 

B = B (z) x + B t (1) 

x z 

( B z a non-zero constant, B x (Oj = 0 ) charged particles trapped in 

the vicinity of the equatorial plane z = 0 where the field intensity is 
weakest do not experience any net magnetic drift in the y direction. In 
the guiding center mode of motion, such particles do experience a non- zero 
curvature and/or gradient drift at any point of their trajectory. However, 
when the effects of such drifts are summed up over one entire bounce period 
between mirror points, it is found that independently of the strength B m 
of the mirror-point field, the net effect of all such drifts cancels out, 
so that the guiding center merely wobbles back and forth in the neighborhood 
of one fixed guiding field line (Figure l) . 


The field (l) is not curl- free and thus is not expected to be formed in 
space plasmas, since a net transport of charge is required for maintaining 
it £ Stem and Palmadessc , 1975] . However, there do exist fields which 
deviate only moderately from (l) -- in the plasma sheet of the earth's 
magnetosphere, in the Jovian magnetodisk, in "spiked helmet" configurations 
in the solar corona ^ Schatten , I971J and possibly in the large-scale structure 
of the interplanetary magnetic field ^ Smith et al. , 1976} • In such "nearly 

drift free" fields the net drii't of particles is not zero but merely small, 
so that a relatively high plasma density can be contained in them without 
creating too high a current density, which might disrupt the magnetic geo- 
metry. Thus nearly drift free fields may be nature's way of containing high-P 
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plasmas effectively, and perhaps this is the reason why they seem to occur 
as frequently as has been observed. 

The purpose of this work lb to analyze guiding center motion in what 
may well be the simplest type of a nearly drift free field, namely the field 

B = | A z x + ( B 0 + 5 k x ) z ( 2 ) 

'■«' d . ~ '-’2 ~ 

where 

k/A = t & 1 

In such a field, as will be seen, .acteristic quantities such as bounce 
times, longitudinal invariants and mean drift rates are readily expressed in 
terms of elliptic integrals, and appendix A develops some of the mathematical 
tools required for handling them. 

One questions which this study answers concerns the dependence of mean 
drift velocities on equatorial pitch angle. For particles in the field ( 2 ) 
confined to the equatorial plane z * 0 the mean drift velocity v^. (which 
equals the gradient drift experienced there — see Stern and Palmadesso [w; 
equation 7 ) is readily computed. How does v vary as the equatorial pitch 
angle is decreased? The results derived here (see Table 1 ) suggest that 

the variation is relatively slow, so that for many purposes results obtained 
for z = 0 can be used as a crude approximation. 

The study also shows that although the average drift velocity of particles 
in such fields is slow, the curvature and gradient velocities observed at 
different points along their trajectory may be quite large. Thus if a distri- 
bution of such particles is observed only in some selected region along their 
trajectories — e.g., near the equatorial plane — misleading values for their 
average bulk velocity could be deduced. 
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THE LONGITUDINAL INTEGRAL 


Consider a somewhat more general two-dimensional field which includes 


both (l) and (2) as special cases 


B ' = B (z) x ♦ B (x) z 

X ' w Z W 


where we assume as before that B x (0) = 0 and that the field is weakest 


at 2 = 0 , so that trapping can occur. Introducing Euler potentials (o< } ji) 


B - V* x 


a convenient choice of these functions is provided by 


oi = / B z dx - / B x dz 


P • 


We now define the longitudinal integral I for a particle with kinetic 


energy W and mirroring field intensity B m by' means of 


J = (2mW) 1//2 § ( 1 - B(z)/B m ) ly/2 ds 


= temW) 1 / 2 I 


with integration performed along a field line between mirror points. Along a 


field line 


and hence 


I = Q>(1- B/B m ) dz 


' 1 


■ 


- 5 - 


However, calculations are considerably simplified if z is replaced as 
integration variable by the field magnitude B . From (j) 

B 2 * B 2 (z) + B 2 (x) (9) 

X z 

Hence along a field line, using (7) and assuming independent variables 
z) 

B dB = [b x dB x /dz + B z (dB z /dx)(-* x/'bz)^ ] dz 

= B x ( d£ x /dz + dB z /dx ) dz (10) 

Thus 

I - f (U) 

J B_ B (dB /dz 4 dB /dx) 

*• A A Z 

THE LIMIT £ 0 

Quantities such as I derived for ( 2 ) will differ only by corrections of 
order £ from their values for the drift-free limit £ * 0 . The calculation 
will thus begin with an analysis of that limit, uting 


B z 

* B o 

= constant 

(12a) 

B 2 

- > 2 

z 2 /k 4 B 2 

(12b) 

B x (dB x /d Z ; 

> - i 

X (B 2 - B 2 ) 

0 

(12c) 

(8/>B q B^ /2 ) 

B m 

f 

(B m - B ) 1//S B 2 dB 

(13) 

J 

B 

0 

( B 2 - B 2 )^ 2 
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fere we have assumed for simplicity Bx(z) * - B x (-z) i in the absence of such 
symmetry, I will be proportional to the sum of two integrals of the form (13), 
one for z > 0 and the other for z < 0 . The calculation then proceeds as 
outlined below and leads to similar conclusions, but I everywhere contains the 
sum of two distinct integrals. In the present case the integral depends on two 
independent parameters B 0 and Bjg , but only one such parameter remains if B 
is replaced by the dimensionless variable u : 


Then 




u - B/B 0 

u o * e iA 

"■ 2 , , 1/2 

(u m ’ u ' 

(u 

1 


(u 2 - 1) 


172 


du 


• (8 b 0 />u; /2 ) s (uj 


(14a) 

(14b) 


(15) 


where is expressible in terms of elliptic integrals K and E of the 

argument m , where m 2 = (u m -l)/(u m 4*l) (see appendix A), namely 


S 2 = (2/15) (ujjj + 1) 1 / 2 [e (2 uj - 9) ♦ K (9 - 2 uj] (A- 10) 


THE HAMILTONIAN W FOR £ = 0 

For some applications it is useful to express the kinetic energy W in 
terms of (j, |i, a procedure which makes it possible to write down 

a Hamiltonian for the twice- averaged motion ^ Northrop and Teller , i960 j 
Chen and Stem , I975 J . Substituting (6) in (15) gives 

j/(2mW) l/2 = (8 B 0 /x4 /2 ) S AuJ (16) 


One now substitutes on the left 


W 


4 B, 


m 


* 


(17) 


giving 


T 


T 


1 


1 
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J n" l/2 - N S^u*) = N S 2 (W/n B 0 ) (18) 

where 

N a 8 (2 m B®)^ 2 /> = constant (19) 

Equation (l8j clearly shove that W Is independent of (o(, f )» a* Is 

required for the drift-free property of the field in the limit 6*0. 

Furthermore, in this limit the dependence of W on n and J may be 

expressed by a function of a single variable -- namely, "W/m depends only 

-l/2 

on the single variable F * J n ' : this property is common to geometries 

free from electric fields, as was shown by Taylor [1964] . 

The functional form of W is defined by the integral relation (6) 
and in order to extract it explicitely that relation must be inverted, 
which tends to be a difficult task. Fortunately, the most significant 
applications of W only require its derivatives, which through the canonical 
equations of motion ^ Northrop and Teller , i960, eqs. 30^ give the bounce 
averaged drift rates and which also, through Liouville 's theorem, allow 
the twice-averaged Vlasov equation to be expressed. These derivatives are 
most easily obtained by implicit differentiation of (6), and this is the 
line which will be adopted in deriving the mean drift velocities when 

It 0 . 

The explicit form of W is required mainly for deriving contours of 

constant Hamiltonian K f Northrop and Teller . 1960^] . In a time- independent 

L orthogonal to B 

magnetic field for which W is known and where an electric field^is present 

and is expressed by its scalar potential , such contours are 

easily drawn and provide a quick visualization of the motion of charged 


f 
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particles through the given configuration. Chen and Stern £l 975 J derived 
the limiting forms of W for a dipole field in the limits J * 0 (equatorial) 
and \i * 0 , and proceeded to bridge these two extreme cases for intermediate 
values of (j, n) with an analytical approximation accurate within 1% . Here 

only the form of W for near-equatorial particles will be developed. 

In the limit J = 0 , (17) gives 

W = 4 B 0 (20) 

Consider next a particle with a small but nonzero value of J ; in that 
case it is useful to introduce a new variable w 


w = U - 1 

(21a) 

= «m- 1 « (W/mB q ) - 1 

(21b) 

S 2 (u m ) = S 2 (w m ) 

(21c) 

sued to denote for convenience 

F = J n -1 / 2 

(22) 

Then (l8) gives 


F/N = S 2 (w/hB 0 - 1) 

(25) 

and one sees that in order to obtain W(n, J) in explicit form it is required 
to invert the function S £ . To lowest order (obtained either by approximating 
the integral in ( 15 ) for small values of v m or by inserting appropriate 
series expansions of E and K in equation A-10) oik finds 

S 2 (w m ) «* i'{2'T]/k) w m 

(Sk) 


r 
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leading to 

W * |iB 0 [l ♦ te'Ji/'jr) f/n] (25) 

Equation (24) also leads to the approximation 

I •* (2 •{? JT /> )(B m - B 0 ) fc6) 

which for the geometry discussed here forms the counterpart of an approxi- 
mation devised by Schulz £1971 > eq. 10 ] for I In the dipole field. 
Using additional terms in the expansions of K and E gives as generali- 
zation of (24) 

s 2 * (/jT JI /4) [w + (7/16) W 2 + (19/256) V 3 ] (27) 

Neglecting the last term and inverting the relation which remains then gives 

W «. n B Q |l + (8/7) [ ( 1 + 1.5756 J u'^/M) 1 ^ 2 - l] j“ (28) 

where N is the constant given by (19). This relation is accurate to about 
1% up to u m *» 1.5 . 
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THE BOUNCE TIME 

The time for one full back- and- forth oscillation, evaluated in the limit 
6, ■ 0 by the same methods as those used for I , is 

fs 

“m 

- (8 B 0 ui /2 /Av) f 

1 

The integral expressed here will be denoted by I' and one easily sees 
that 

I' = 2 dS 2 /du m (30) 

(compare Schulz W , eq. 7 h To evaluate I* it is best to 
split it into two parte 



• I i + I 2 (31) 

Integrating by parts, one finds that I x reduces to one of the class of 
integrals analyzed in appendix A , namely 

I x - 2 S 1 (32) 

The second integral may be handled by standard substitutions whi :h lead 
to elliptic Integrals (e.g. Milne-Thomson [las*] , page 597 ) and is 
given by 

"" l ' T !ll! - ‘ 



I 
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l/p 

I a ■ 2 (% *l) K(m) (33a) 

where K is a complete elliptic integral of the first Rind and its 
argument ra (which also appears in the termB of eq. 32) satisfies 

m 2 * (u^ 1 )/ 1 ) (33b) 

l/p 

A tabulation of 1=2^ I' is given in table i . The conversion of 
these values to actual time intervals requires multiplication by 4 B Q />v 
in conditions appropriate to the earth's magnetotail 

E c /> * 1 R e 

and for 1 keV protons T turns out to be approximately equal to T in 
minutes. For 1 keV electrons, 1.36 T approximates X in seconds. 

NORMALIZED CALCULATION FOR tf- 0 

Before calculating the motion of charged particles in the field (2) , 
it is useful to transform the equations to dimensionless form, since that 
reduces the number of independent variables and therefore allows the 
same result to be used for a large variety of situation . Such transfor- 
mations are usually accomplished by scaling -- e.g. dividing all distances 
by a scale length. In the present system a natural scale length is 

L * B 0 /A (34) 

One can visualize L a6 follows. In the limit £,=0 ary field line 
follows a parabola 

o( ■ Az 2 /4 B 0 x = constant (35) 
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Along this parabola the direction of B gradually rotates from being 
orthogonal to the x axis to being parallel to it. The point at which B 
has rotated by 45° — half its total rotation range -- occurs when B x = B z 
and is L/2 above or below the equatorial plane. As noted at the end of 
the preceding section, in the earth's magnetotail L *# 1 Rg . We now use 


L to define scaled distances 

^ = s/L (36a) 

5 = x/L (36b) 

and also introduce a scaled magnetic field 

u - B / B 0 (37a) 

and a scaled Euler potential 

o<' = */B 0 L (37b) 

The dependence of ^ and w on ^ is depicted in Figure 2 . 
With these variables, the basic relations become 

(58) 

ot' - 5 ♦ e$ 2 A - S 2 A (39) 


^3B x / < 4) z = X ^u x /^^ 

J'lx - > Du z /T>$ 

Substituting in (ll ) 

I „ kL f (1 ' Am) 1 ^ 2 u a du 

~ u z u x (du x /d$ + du z /d$ ) 
u(0) 


(40a) 

AOb) 

(U) 


T 


r 




Prom ( 38 ) the sum in the denominator is simply (l + £)/2 . The components 
of £ are also easily derived, but their functional dependence must be trans- 
formed from ( tj) to (u, ^') in order to suit the integration. If terms of 
order £ 2 or smaller are neglected, one finds 

u a * Tj 2 /4 + 1 ♦ 

« 'f/k fl+E( + % z /k) (42) 

(l + £) u£ » u 2 - 1 - 

u x « (u a - 1 - t*') 1 / 2 (l + E.)' 1 / 2 (43) 

From (38), (39) and (43), using ( 39 ) only to lowest order 

u 2 * 1 + e 5/2 

* 1 + ( E/2)( <X’ + u|) 

* 1 + ( + u 2 - 1) (44) 


Substituting everything in (4l) one obtains the form of I to order £ as 
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Note that the dependence on must he preserved even If we choose to 


deal with the 

field line oi* “ 

0 , since n)l/ r k<ai t 

is required in the 

calculation of mean drift velocities. Introducing a new variable w (not 

related to w 

in 

equation 21 ) 

through 




w* . 

u 2 - 

(4?a) 



V 2 * 

m 

4- &* 4 

(47h) 

we get 







u 

w [l +£(o< , / 2 v 2 )] 

(48) 



U m - u « 

(w m - w) [jL - £( ot'/2 

v w m )J ( 49 ) 

Thus to order € 




ii 

M 

f 

(v m . u) 1 / 2 V 2 

f 1 + £( ^ ‘ 

L 2 w 2 

o( oi 4- w 2 - 1 

J 

1 

(» 2 - 1) 1/2 

4 w W m 2 




dw 


= ~ 4 j 2 [ S 2 (w m> + £(<*' S 0 /2 - *(' Sjk v m - (*T - 1) S s /2 - S^/2 ) J 


where the are Integrals defined in appendix A 

“n / v 1/2 

r I Mm- U \ 

Sk(%) 


• r fer) 


k 

u du 


Now let I ' he defined by 
I 


- 1/2 TI 
u w 1 


(50) 


(A-l) 


(51) 


iilii . MM 
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By (6) 


where by (47b) 


J * (2 n |i Bq) 1 / 2 I'CaC', v m ) 


w m “ w m^ u m* d') = w m (wA© 0 > «< ' ) 


The last two relations implicitely express the functional relationship 


W - w(o< , n, J) 


Differentiating by 1 


HU' ^1' 

'Dot n>e( 


i + 7)1' / 1 Ow 

“ \ HB 0 7) u 


Only errors of higher order are committed if we set both ^ and 'Dw E j/^u ffl 
equal to 1 and if we evaluate (in the expression below) *7)r'/^w only 
to ihe lowest order. This gives 


5 lS J du m 


2r . dS 2 ^ w m 

7)o(' du m 


From (48) 


hence 


*7) " t/2 u m 


£ u L 


j^S 0 /2 - Sjh u m - S 2 /2 - ( dS 2 /dUffl)/2 UnJ ( 58 ) 
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SCALING OF BRUT VELOCITY AND TIME 

An independent numerical integration of guiding center drifts in the field 

(2) was performed which served as a check 

which 

on the analytical formula (58) and^also provided information shout the 
drift velocity in various portions of the orbit and about the time spent 
by the particle in these portions. In performing these numerical checks some 
additional scaling was introduced, in the following way. 

The mean drift velocity <Vy> may be derived from (58) by the relation 

« <V> ' - • H? f(u m ) <59) 

where f(u m ) is the function defined by (58) ; note that since 'DW/U** is 
expressed in terms of the derivatives of I (see eq. 55) , we are actually 
using equations (27) of Northrop and Teller [i960] , rather than equations 
(50) there. Let a scaling velocity v Q be defined 

v Q = 4km v^q (60) 

As will be seen below (equation 72), v Q is the gradient drift velocity for 
equatorial particles (j = 0 ), also given by equation (7) of Stern and 
Palmadesso £l975j • By (59) 

= - 2 v 0 f(u m )/u ffl (61) 

One may also devise a "scale time" t Q for the bounce motion 

t 0 = L/v = B 0 />v (62) 


so that (29) becomes 
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% 

T - 8 t c 4 /2 j u s [( V uKu 2 - 1)] ' 1/2 du (63) 

1 

Table 1 lists < v y ^> /v q and T/4t 0 for selected values of u m , obtained 
by means of ( 58 ) and (63), and Figures 3 and k represent these functions 
graphically. 

Let v (vithout averaging brackets) be the guiding center drift velocity 

v 

of particles with given v Q and , computed for various points of the 

trajectory. For such particles 

<v y > r - §Vy at (64) 

This relation allows <v > to be derived numerically, but one must proceed 
with caution, because the integral has to be evaluated to an accuracy of the 
first order in £ . Transforming the integration variable by means of a 
zero-order relation between t and seme other variable (e.g. u or z ) 
may lead to a finite error even if v in the integrand is given correctly 
to order £ . One may, however, us ( 7 ), which is exact, to obtain 

<\> - 


where z ffl corresponds to the mirror point. Note than at the i^per limit 
the integrand diverges, so that a different mode of numerical integration may 
be required for the vicinity of the mirror point: one convenient device for 
accomplishing this is described in appendix B . 
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THE DRIFT VELOCITY 


All guiding center driftB in the field (2) are in the y direction and 

A 

therefore can he added or subtracted algebraically. This follows from the 
observation that the components constituting these drifts (see below) are 
aligned either with B * VB 2 or with V * B . In the first case 


B * V B* 


y ( - X 3 z 2 A + B? k ) 


= y aBq ( • $ 2 A e ) 


and in the second one 


V x B = y ( h/2 ) (l - £ ) 

w 

The two principal components of v w are the gradient drift 

+* y 

V = (v 2 /2u>B a )(B X VB) 

A»g * #«* 

and the curvature drift 

v c = (v 2 /u) B) Bs(B-VB) 


= (v 2 /cJB) j^(Vx B) + B X 7bJ 


where U) = <£B/m is the gyration frequency. Then by (60) 


v^ /2u)B 2 = 2 v Q /(k % u B) 


v 2 /(JB 


4(v 0 /ku 2 )(l - u/ujj) 


v o 1 / t 2 , f X 

v g * 7" ' ‘ "IT * ^ ' 

6 e. u a u 3 4 


1 


1 
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These results may now be Inserted In ( 65 ) to obtain a numerical estimate 
of <v y > . Note that although v y may be separated in ( 72 ) and (73) Into 
parts of order v Q and v 0 /g. , the contribution of the latter part to 

^v y > does not vanish but is merely of order v Q , because the integration 
path differs slightly from its limit for £= 0 . This prevents one from 
calculating through ( 65 ) by summing up terms of the same order and 

instead, one must include terms of lower order , the contributions of which 
almost (but not completely) cancel out. 

This in its turn demands a computing accuracy exceeding by a factor l/fi 
the accuracy of the final result: the advantage is that a relatively simple 
formula ( 65 ) can be used, avoiding the cumbersome transformations which head 
to ( 58 ). In practical evaluations, the numerical result agreed within a few 
percent with ( 58 ); as an additional ohecls, the time was also integrated 
numerically, and the results were compared with ( 29 ) • 

Table 2 contains a sample result for a field with fc,= 0.01, for particles 
mirroring at 1.4 , u m = 1.2227. The fourth column tabulates values of 

the integrand of ( 65 ), multiplied by the step size Az = 0.05 (for brevity, 
only every alternate step is tabulated), while the 5 th column gives the inte- 
grand of (63 ) , multiplied by 2 <Az . Analytical and numerical results 

for the total displacement and time (equator to mirror point) are 
listed, and the displacements obtained analytically are shown in Figure 5 • 

Note that the displacement is negative near the equator, where curvature drift 

dominates, but is positive near the mirror points (the net drift is also posi- 
tive). Note also that the drift’ velocity deduced for the point at the top of 
the table is about 50 times larger, and in the opposite direction, than the 
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A certain amount of caution should be excercised when the preceding results 
are applied to observations, because the average drift velocity will not in 
general equal the plasm's bulk velocity orthogonal tc the magnetic 

field B . To emphasize this point, consider the current density 1 carried 
in a direction orthogonal to B by a plasma of particles of a single kind 
only, e.g. electrons. If at any given point f is the distribution f.inction of 
such particles and n their density, then their bulk velocity is 

Vj, = (l/n) fv 4 f d 3 ^ (74) 

r* J <v 

and 

j = n q V* (75) 

On the other hand, if is the mean guiding center velocity and M 

the magnetic moment density, then 

j = n q <v.> + V X M (76) 

Unless V#M vanishes <>£> will in general differ from V± .In particular, 
one can envision static equilibria in which the plasma distribution function is 
everywhere isotropic: by (74) V A and both vanish, but <v^> may be 
quite large. 


As a concrete example, let B - B(x, y)z,n = constant and let all 
particles have the same v^ . Then no curvature drift exists and all particles 
share the same gradient drift velocity, given by (68). However 


M 

V * .¥ 


n m vf 
2 B 

n m v 3 
2 B 2 


A 


B 


S X 7B 


(77) 

(78) 


and J. of (74) vanishes (this still holds if populations of this kind with 

M * 

v^ are superposed). Thus V = 0 while clearly <v^> / 0 . 


different 
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THE OBSERVED VALUE OF £ 


L -. from 

19741 have concluded observations that at 30 Rg 

the average thickness of the plasma sheet is *» 2.5 Rg # over which B x changes 

from 10 7 to - 10 7 . This implies 


^ = 2 *» 16 7/ Rg 


(79 ^ 


On the other hand, Behannon [woj reports that at 30 Rg , B z *• 1.87 7 » 
and that the average value of B z changes from x = -20 Rg to x = - 40 Rg 
(x being negative on the nightside) by about 2 .5 7 j giving 


k = 2 7)B 2 /'Dx «** 0.25 7/Rg 


(80) 


Thus t *»» l/64 , which means that the order of approximation used here is 
indeed an appropriate one for the geomagnetic tail at 30 Rg • At x » -20 Rg 
the plasma sheet is about twice as thick, suggesting X is close to half 
its above value, while k ^ Behannon , 1970^J is about twice the value 
of ( 80) : £ is thus considerably larger, but the approximation will still hold. 


GENERALIZATION 

One can in principle extend the preceding calculation to more general 
fields of the nearly drift-free type, e.g. to fields with as in (2) 
but with B x expanded in a more general fashion 


B x (z) 


o > 2 * c ; 


z 2 + c. 


( 81 ) 
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where the series is presumed to converge and where even-numbered coefficients 
vanish in a symmetric configuration. The complexity of the calculation, 
however, increases rapidly with the addition of terms, and the use of 
numerical formulas such as ( 65 ), with appropriate generalizations of ( 72 ) 
and (73 )> may well be quicker and less prone to errors than an analytical 
approach. Here we shall sketch out one relatively simple generalization — the 
derivation of the bounce time X* in drift-free fields (B 2 * B Q = const.) 
where B x has the more general expansion (8l). To keep the calculation short 
we will also assume that all even powers in (8l) vanish, so that all bounce-time 
integrals reduce to 4 times the integral between z = 0 and the mirror point 


on either side. Thus 


+ B o! 


Dividing by B Q and introducing new variables 


* i 


By ( 7 ), and with the notation used earlier, this gives 


T - 


ds/v M 


= (4 u^ 2 / v ) \ (ujjj- u) 1//2 u ('tJz/'hu)] du ( ) 


Let the function B x (z) in (8l) have an inverse 


z(B x ) 


Z. a * B > 


where again in the symmetric case k only takes odd values, since z and B x 
reverse signs together, and where = 2/> . Defining L and £ as in (34), ( 36 ) 


1 
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■ L i > k/2 

(85) 

where 



** = 


(86) 

Thus 



dz/du 

- (l/>) 4 $/ 4 u “ £ k u (u 2 - l) lk/,Z 

“ lJ (8?) 

Defining t as 

in (62) 


T 

■ kt o“i /2 72 * \ w 

k=l 

(88) 

where 

T k 

% 

= ^ (u m - u)’ l/2 u 2 (u 2 - i } (^“2 ) /2 du 

(8S) 


1 


Because k is odd, any (except for T , which appears 

in (29) and 

which has been derived earlier) may be expressed as a linear combination of 

integrals of the type 


v k 

“m 

= ^ [(u 2 - l)/(u m - u) J U k du 

( ) 


1 


© *8 • 

T 3 

* 


T s 

. V - V 
4 S 


By equations (A-8) and (A-6) these integrals can be expressed in terms nt 


the integrals of appendix A , using the relation 
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V k * [(^ + 5/2 ) + 2 ' (k ♦ 3/2) (91) 

If even-numbered values of k are also allowed In the expansion of B x , 

the calculation * s modified somewhat and integrals of the type (89) 

corresponding to even values of k also appear. Such integrals do not have 

to be expressed in terms of elliptic integrals but instead may be integrated 

in closed form after a substitution v = u m - u , since their integrand is 

- 1/2 

then a polynomial in w multiplied by w ' 

CONCLUSION 

This work derives for the first time the properties of adiabatic motion in 
a tail- like two-dimensional magnetic configuration possessing a finite 
/ 'lB z /'' 5 x , which is essential for a net cross-tail drift. It derives 
instantaneous and average drift velocities, bounce times, longitudinal invariants 
and approximations to the adiabatic Hamiltonian, and the results are suitably 
scaled to reduce as much as possible the number of independent parameters. 

The mean drift velocity is (as expected) much smaller than the instantaneous 
drift velocities encountered at typical points on the trajectory. It ic com- 
parable to the drift velocity v^ of particles confined to the equatorial 
plane, which are sensitive only to the small gradient c m B z : as pitch angles 

decrease, this velocity becomes smaller and for the range investigated here it 

pairs of 

levels off around v Q /j • It 1 b interesting to note that if^ neutral points and 
strong dawn-to-dusk electric fields are formed in the tail during cub storms (as 
is widely believed), the gradient of B z will reverse its direction on the 
anti-sunward side of such points and particles drifting there will lose energy 
to the electric field. This, however, is only expected in r. very limited region. 

! - ! ] 1 t — 


I 


I 


1 


I 
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This slow drift suggests that In the plasma sheet of the magnetospherlc 
tail particles advance in the dawn-dusk direction rather slowly, at speeds 
of the order of 1-2 Rg/hour . If an electric field of, say, 40,000 volts 
exists across the sheet, the earthward convection of such particles will he 
considerably more rapid. The present calculation does not apply to non-adiabatic 
protons, which constitute an important part of the plasma sheet's population, 
but because the present geometry differs only slightly from one which is 
drift-free even for non-adiabatic orbits £ Stern and Palisades so ^l^il , the 
cross-tail motion even then is expected to be relatively slow. 

It should be readily possible to extend this work to the magnetodisk 
configuration observed near- Jupiter. More generally, the results derived here 
make available for a tail- like configuration quantities which are already 
well-known for the dipole case and should make possible the generalization 
of calculations for the dipole geometry to such configurations sb well. 


Appendix A : A CLASS OF BfTEGRAIS 


In the study of particle motion in newly drift free configurations one 
frequently encounters integrals of the fom 

W - 

These integrals and their derivatives are readily expressed in terms of 
complete elliptic integrals K(m) and E(m) of the first and second kind 
£e.g. Milne - Thomson , 1964 J , of the argument 

m = r^ u m“ 1 ^'V’ ^ ^ (A-2 ) 

In what follows , unless otherwise is stated, K and E will be understood 
to be evaluated for the above argument. The first two integrals of this type, 
corresponding to k = 0 a: 1 , may he expressed by standard substitutions 

or by consulting tables ( Gradshteyn and Ryzhik £ 1965 J , integrals 3.141.4, 
3.141.22 and 3-141.23); in either case cue finds 

S Q * 2 (u^+1) 1 / 2 (K - iZ) (A- 3) 

S x = (2/3) (u^l) 1 / 2 ( u m E - K ) (A- 4) 

For higher values of k a recursion formula may be developed. Let 

\ = / n* 2- 1)( v x/ " 2 uK dj (A-5) 

where the implied limits of integration, here and in subsequent formulas, are 
always those of (A-l) . Multiplying numerator and denominator by (u“- l) then gives 
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Sc -2 S k S k -2 ^ A-6 

On the other hand, integrating (A-l) by parts 

s k - J ( % - U) 1 ^ 4 [( u *. D^j 

- - <* • 5> v 2 + (u m/2) Va (A ’ 7) 

where is defined in ( 90 ) and may be integrated by parts 

\.g * * 2 J <“*- 1) » k ' 2 d[(v u)l/2 ] 

= 2 <* * 21 «fc-5 + 2 Vl (a-8) 


Substituting the last three equations in a way that eliminates all Integrals 
not of the form (A-l) yields the recursion 

i **! 1 S * ‘ - 1) s k -l *(3L-\) S k . 2 - u D (k - 2) S k . 3 (A-9) 

For k = 2 the last term drops and one obtains 

S 2 = (2/5) u m S x + (3/2) S 0 

= (2/15) (i^tl) 1 ^ (2 u| - 9 ) +• K (9 - 2 uj] (A- 10 ) 

To obtain dS 2 / du m one differentiates beneath the integral sign, noting 
that the contribution from the upper limit vanishes. The result is equation 
( 30 ) which has already been analyzed, giving 

dS 2 / du m - S 1 + (u m + l )' 1 / 2 K 

1/2 % E + 3(u m -»-l)' :L/2 (l - 2 Uj,,) K (A- ID 


(2/3) (u m + l) 
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Performing the same steps for arbitrary k 

d^/d^ = 1/2 J u k [(u 2 - l)(u m - u)]’ l/2 du 


“ ( 1 / 2 ) 


(A- 12) 


One can now express V. by means of (91) and thus obtain a recursion 
for dS^/dUjjj . For even values of k the recursion begins with (A-il)j 
for odd values one needs dS J du m , which can be obtained either by 
standard substitutions or from tabulations ( Gradshteyn and Ry^hik £1965} 
integral 3. 132. 5 ). In either case 


dSj/dUjn - (u m + l) 1 / 2 E - ( % + I)' 172 K 


(A- 13 ) 
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Appendix B : INTEGRATION REAR u = u ffl 

At the mirror point v # = 0 and therefore special care is needed when (29) 
or (65) are integrated near that point. The problem resembles the evaluation 
of integrals of the form 

A 

f f(x) x' 1 / 2 dX (B-l) 

0 

near their lower limit and there exist various advanced methods for handling 
it ^ Smith and Bewtra, 1976 J . In this work a rather crude but effective 
approach is used, as follows. 

Consider the evaluation of (B-l): one first divides the range into two 
parts, the main range B X A which Is handled by conventional methods 
(e.g. Simpson's rule) and the end region 0 ^ X < B . In that region f(x) 
is approximated by a polynomial 

f(x) * f (0) + a X +■ b X 2 (B-2) 

from which the region's contribution to the integral is found at once to be 
very nearly 

2 B^ 2 ( f(0) taB/3 + b B 2 /5 ) (B-j) 

The coefficients a and b can be evaluated using f(B) and some intermediate 
value f(C), 04C<B , but computing time is saved if one directly evaluates 

aB and bB 2 , which appear in (B-3 ) • An alternative approach first removes the 

l/2 

singularity in (B-l) by transforming to a new variable w = X . If f is 
then expanded in even powers of w , (B-3) is obtained; if odd powers are also 
allowed, the result is the same as when both integral and half integral powers 
are used in (B-2). 


The integrals ( 29 ) and ( 65 ) readily reduce to the fora (B-l) if one sets 
X = u m - u . Note that over the end region (hut not over the min range) ( 65 ) 
must first be transformed into an integral over u , after vhich one introduces 
dX = - du ; by (42) and (36a) this requires the addition of a factor 



(B-4) 


The above factor contributes to f(x), but since f(x) is used only 


numerically, z does not have to be expressed in terms of u and the above 
factor can be used in the form given in (B-4). 



■ 31 - 


REFERENCES 


Behannon, K.W., Geometry of the Geomagnetic tail, J. Geophys ■ Res . , Jp , 

743-753, 1970. 

Bovling, S.B. and R.A. Wolf, The motion and magnetic structure of the plasma 
sheet near 30 Rg, Planet . Space Gel . , 22 , 673-686, 1974. 

Chen, A.J. and D.P. Stern, Adiabatic Hamiltonian of charged particle motion 
in a dipole field, J. Geophys . Res . ,80, 690-693, 1975* 

Gradshteyn, I.S. and I.M. Ryzhik, Table of Integrals , Series and Products , 
translated from the 4th Russian edition by A Jeffrey, Academic Press, 

1965- 

Milne -Thomson, L.M., Elliptic integrals, in Handbook of Mathematical Functions , 
M. Abraraowitz and I. A. Stegun, editors, U.S. Govt. Printing Office, 

1964, pp. 587-626. 

Northrop, T.G. and E. Teller, Stability of the adiabatic motion of charged 
particles in the earth's field, Phys . Rev . , 117, 210-225, i960. 

Schatten, K.H., Current sheet magnetic model for the solar corona, Cosmic 
Electrodynamics , 2, 232-245, 1971* 

Schulz, M. Approximate second invariant for a dipole field, J. Geophys . Res., 

J6, 3144-3148, 1971. 



* 32 " 

Smith, E.J., B.T. Tsurutani and R.L. Rosenberg, Pioneer 11 cfeaerrations of 
the interplanetary sector structure up to 16° heliographic latitude 
(abstract SS-8), EOS, £2* 997, 1976. 

Smith, P.H. and N.K. Bevtra, Dependence of charge exchange lifetimes on 
mirror latitude, Geophys . Res . Lett . , 689 - 692 , 1976. 

Stern, D.P- and P. Palisades so, Drift- free magnetic geometries in adiabatic 
motion, J. Geophys. Res ., 80, 4244-4248, 1975 • 

Taylor, J.B., Equilibrium and stability of plasma in arbitrary mirror fields 
Phys . Fluids , J, 767-773, 1964. 


- 35 - 


CAPTIONS TO FIGURES 


Figure 1 -- Schematic view of the motion of a charged particle in a drift-free 
magnetic field of the type discussed in this work. 

Figure 2 -- The configuration of a magnetic field line (left scale, lower 
curve) and the variation along it of u = B/B 0 (right scale, 
upper curve) in the limit £,= 0 . Both x and z are normalized 
in accordance with (56). 

Figure 3 — The variation of the mean drift velocity, in units of the equa- 
torial drift velocity v Q , as a function of u = b/b 0 . 

Figure 4 — The quarter-hounce time, in the normalized units of Table 1, 
as a function of u = b/b o . 

Figure 5 — The field line of Table 2, divided into segments each of which 

corresponds to A z = 0.05 (normalized units). Drawn at the end of 
each segment is a line proportional to the net distance covered 
by guiding center drifts during the particle 's traversal of the 
segment (actually, the displacement should be orthogonal to the 
drawing). The net displacement for the entire quarter period is 
about equal to the lowest displacement drawn (but in the opposite 
direction); the time spent traversing that first segment, however, 
is only about I/ 5 O of the quarter- period. 
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